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The ferromagnetic phase of an Ising model in d = 3, with any amount of quenched antiferro- 
magnetic bond randomness, is shown to undergo a transition to a spin-glass phase under sufficient 
quenched bond dilution. This general result, demonstrated here with the numerically exact global 
renormalization-group solution of a d = 3 hierarchical lattice, is expected to hold true generally, 
for the cubic lattice and for quenched site dilution. Conversely, in the ferromagnetic-spinglass- 
antiferromagnetic phase diagram, the spin-glass phase expands under quenched dilution at the 
expense of the ferromagnetic and antiferromagnetic phases. In the ferro-spinglass phase transition 
induced by quenched dilution reentrance as a function of temperature is seen, as previously found 
for the ferro-spinglass transition induced by increasing the antiferromagnetic bond concentration. 

PACS numbers: 75.10.Nr, 64.60. ah, 75.50.Lk, 05.10.Cc 



The spin-glass phase [l[ is much studied due to its 
prominent role in complex systems, as an example of 
random order. In its simplest realization in the Ising 
model, the underlying system has randomly distributed 
ferromagnetic and antiferromagnetic bonds. In spatial 
dimension d = 3, at low temperatures, ferromagnetic or 
antiferromagnetic ordered phases occur when the system 
has predominantly {e.g., more than 77% [^) ferromag- 
netic or antiferromagnetic bonds, respectively. In be- 
tween, the spin-glass phase occurs. The occurrence of 
the spin-glass phase, in which the local degrees of freedom 
are frozen in random directions, has strong implications 
in physical systems that are realizations of the spin-glass 
system, spanning from materials science to information 
theory and neural networks. 

We have studied possibly the simplest modification of 
the spin-glass system, to be commonly expected or real- 
ized in physical systems, namely, the removal of bonds. 
We find important qualitative and quantitative effects. 
This Ising spin-glass system with quenched bond vacan- 
cies [3,,,^^^, J, 8, 9, 10] has the Hamiltonian 

-/3i/ = (1) 

where s; — ±1 at each site i and (ij) indicates sum- 
mation over nearest-neighbor pairs of sites. The lo- 
cal bond strengths are Kij — K > with probabil- 
ity p+, Kij — with probability or Kij — 
with probability q = 1 — — p_ , respectively corre- 
sponding to a ferromagnetic interaction, an antiferro- 
magnetic interaction, or a bond vacancy. This model 
has previously been studied at zero temperature [1, Q 
and in its spin-glass phase diagram cross-section [l^l by 
position-space renormalization-group theory, in n-replica 
version at zero-temperature 0,1a], by mean-field theory 
[3] , and by momentum-space renormalization-group the- 
ory around rf = 6 dimensions [i^] , and by series expansion 
fc^]. 

We have performed the numerically exact 
renormalization-group solution of this system on a 



d = 3 hierarchical lattice [HI, [13, [13] > be given below. 
Exact solutions on hierarchical lattices constitute very 
good approximate solutions for physical lattices. |14| 
We calculate the global phase diagram in the variables 
of temperature l/K, bond vacancy concentration q, 
and antiferromagnetic bond fraction p-l{jpj^ -I- 
obtaining a rich structure and finding reentrance as 
a function of temperature induced by bond vacancy. 
Our results agree with and extend the previous results 
[3, 4, 5, 6, 7, 8, 9, 10] on this system. 

Our results are most strikingly seen in FiglT] The top 
curve in Fig. 1(a) corresponds to the quenched dilution of 
the system with no antiferromagnetic bonds (p_ =0). As 
the system is quench-diluted, by increasing the missing- 
bond concentration q, the transition temperature to the 
ferromagnetic phase is lowered from its value with no 
missing bonds at g = 0, until it reaches zero temperature 
and the ferromagnetic phase disappears at the percola- 
tion threshold of g = 0.789 (to be com par ed with the 
value of 0.753 in the simple cubic lattice flo]). However, 
with the inclusion of even the smallest amount of anti- 
ferromagnetic bonds (lower curves), a spin-glass phase, 
extending to finite temperatures, always appears before 
percolation. This result was previously obtained at zero 
temperature [1, H, 0, 9 and around d = 6 [4] . 

The phase boundary between this vacancy-induced 
ferromagnetic-spinglass phase transition shows reen- 
trance, as also seen [l^, [13, [Ilj in conventional spin- 
glass phase diagrams where the antiferromagnetic bond 
concentration is scanned. In Fig. 1(b), where the curves 
correspond to higher percentages of antiferromagnetic 
interactions among the bonds present, starting with 
P-/{p+ + P-) — 0.25 in the top curve, the ferromag- 
netic phase has disappeared and only spin-glass ordering 
occurs. As seen from Fig. 1(a), the percolation thresh- 
old of the spin-glass phase is slightly lower than that of 
the pure ferromagnetic phase and, before the disappear- 
ance of the ferromagnetic phase, the percolation thresh- 
old of the spin-glass phase has a slight dependence on 
p_/ (p++p_). The percolation threshold of the spin-glass 
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FIG. 1: (Color on-line) Calculated phase diagrams at constant 
p_/(p+ +P-). All phase transitions (full lines) are second or- 
der. The dashed horizontal lines (shown only partly in (a)) 
are not phase boundaries, but the Nishimori symmetry lines, 
given by Eq.((2]|, for each value of p-/{p+ + p-). The multi- 
critical points in (a), mediating the ferromagnetic, spin-glass, 
and paramagnetic phases, occur on the Nishimori symmetry 
lines. These lines do not cross the phase boundaries at any 
other type of point. Thus, the Nishimori symmetry lines in 
(b) are at temperatures above the phase boundaries. 

phase settles to the value of 0.763 after the disappearance 
of the ferromagnetic phase. 

Fig. 2 shows the conventional phase diagrams of tem- 
perature versus the fraction + p-) of antiferro- 
magnetic bonds in the non- missing bonds, at fixed values 
of the dilution q. As the dilution is increased, the phases 
are depressed in temperature, as can be expected. How- 
ever, simultaneously, it is seen that the spin-glass phase 
expands along the P-/{p+ axis, at the expense 
of the ferromagnetic and antiferromagnetic phases, even- 
tually dominating the entire low-temperature region. It 
is also seen, for q = 0.763 and higher, that the spin-glass 
phase occurs in the phase diagrams as two disconnected 
regions, near the ferromagnetic and antiferromagnetic 
phases. A similar disconnected topology has recently 
been seen in the Blume-Emery-Grifhths spin glass. [l9] 

The dashed lines in Figs.l and 2 are the Nishimori 
symmetry lines [13, [HI , 

e±2^ = ^ . (2) 
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FIG. 2: (Color on-line) Calculated phase diagrams at con- 
stant quenched dilution q = 1— p+ +p-. The phase diagrams 
being symmetric about p-/{p+ + p-) = 0.5, with the anti- 
ferromagnetic phase replacing the ferromagnetic phase, only 
the p-/{p+ + P-) < 0.5 halves are show. All phase transi- 
tions (full lines) are second order. The dashed curve is the 
Nishimori symmetry line given by Eq.((2j. All multicritical 
points, mediating the ferromagnetic, spin-glass, and param- 
agnetic phases, lie on the Nishimori symmetry line. 



led system are on the Nishimori symmetry lines [22l . |23| , 
as also previously seen [13] for this system. Thus, as 
illustrated in Fig. 2, it is possible to continuously popu- 
late, with multicritical points, the low-temperature seg- 
ment of the Nishimori line, by gradually changing the 
quenched dilution q. The Nishimori symmetry condition 
appears in Fig.l as a horizontal line for each value of 
P-/{p+ +p_). This horizontal line intersects the up- 
per curve in Fig. 1(a) at zero temperature, thereby im- 
plying the occurrence of a zero-temperature multicritical 
point at the percolation threshold, as also deduced from 
the sequence of phase diagrams in Fig. 2. The horizontal 
lines of the Nishimori condition intersect the two other 
phase diagrams in Fig. 1(a) at their multicritical point. In 
Fig. 1(b), multicritical points do not occur and the hor- 
izontal lines of the Nishimori condition do not intersect 
the phase boundaries, occurring at higher temperatures 
than the phase boundaries. 

Fig. 3 shows the zero temperature limit of the global 
phase diagram of the currently studied system. In the 
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FIG. 3: (Color on-line) Calculated zero-temperature, 1/K 
0, phase diagram. All phase transitions are second order. 




hierarchical lattice of Fig. 4 has been used in this work, 
because it gives numerically accurate results for the crit- 
ical temperatures of the d = 3 isotropic and anisotropic 
Ising models on the cubic lattice. [14] 

In systems with quenched randomness, the 
renormahzation-group transformation determines 
the mapping of the quenched probability distribution 
V{K)\?,S\ At each step, the innermost unit of the lattice 
as pictured on the right side of Fig. 4 is replaced by 
a single bond. This is effected by a series of pairwise 
convolutions of the quenched probability distributions, 

V{k)^ J dK^dK^^Vi{K^)rii{K")S{K~^R{K^,K")), 

(3) 

where R{K^,K'') is 



K. 



K, 



II 



(4) 



for replacing two in-parallel random bonds with distribu- 
tions Pi{K') and Vii^K^^) by a single bond with V{K), 
or 



(5) 



FIG. 4: The d = 3 hierarchical lattice for which our calcula- 
tion is exact is constructed by the repeated imbedding of the 
graph as shown in this figure. This hierarchical lattice gives 
numerically accurate results for the critical temperatures of 
the d = 3 isotropic and anisotropic Ising models on the cubic 
lattice 

zero-temperature phase diagram, it is again seen that 
a spin-glass phase intervenes d, H, 0, @], with smallest 
amount of quenched antiferromagnetic bonds, between 
the ferromagnetic phase and percolation, causing a direct 
ferromagnetic-spinglass phase transition. 

Our method, detailed in other works [H, [13, [H, [ll] , 
will be briefly described now. We use the d = 3 hi- 
erarchical lattice whose construction is given in Fig. 4. 
This hierarchical lattice has the odd rescaling factor of 
6 = 3, for the a priori equivalent treatment of ferro- 
magnetism and antiferromagnetism, necessary for spin- 
glass problems. Hierarchical lattices admit exact solu- 
tions, by a renormahzation-group transformation that re- 
verses the construction steps. [Ill, Thus, hierarchi- 
cal lattices have become the testing grounds for a large 
variety of coo per ative phenomena, as also seen in recent 

works, dl MM, [sO, ill, [13, 111 H HE m, 113 . The 



for replacing two in-series random bonds by a single 
bond. The probability distributions are in the form of 
probabilities assigned to interaction values, namely his- 
tograms. Starting with the three histograms described 
after Eq.([T]), the number of histograms quickly increases 
under the convolutions described above. At a computa- 
tional limit, a binning procedure is used before each con- 
volution to combine nearby histograms [H, [H, [H, HI] , 
so that 160,000 histograms are kept to describe the prob- 
ability distributions. The flows of these probability dis- 
tributions, under iterated renormalization-group trans- 
formations, determine the global phase diagram of the 
system. 
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